Precise MS light-quark masses from lattice QCD in the RI/SMOM scheme 
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We compute the conversion factors needed to obtain the MS and RGI up, down, and strange- 
quark masses at next-to-next-to-leading order from the corresponding parameters renormalized in 
the recently proposed RI/SMOM and RI/SM0M 7jJ renormalization schemes. This is important for 
obtaining the MS masses with the best possible precision from numerical lattice-QCD simulations, 
because the customary RI^/MOM scheme is afflicted with large irreducible uncertainties both on 
the lattice and in perturbation theory. We find that the smallness of the known one-loop matching 
coefficients is accompanied by even smaller two-loop contributions. From a study of residual scale 
dependences, we estimate the resulting perturbative uncertainty on the light-quark masses to be 
about 2% in the RI/SMOM scheme and about 3% in the RI/SM0M 7(J scheme. Our conversion 
factors are given in fully analytic form, for general covariant gauge and renormalization point. We 
provide expressions for the associated anomalous dimensions. 



Lattice QCD has, in recent years, seen important 
progress on several fronts: there exist lattice regulariza- 
tions preserving exact chiral symmetry in the limit of 
vanishing quark masses, while algorithmic and techno- 
logical advances have put lattices fine enough to simu- 
late physical light-quark masses within reach. As a re- 
sult, nonperturbative results in the physics of light quarks 
with a precision of a few percent or better become achiev- 
able with current or upcoming simulations [1] . These in- 
clude the masses of the light quarks, as well as hadronic 
matrix elements such as B^, figuring prominently in 
the unitarity-triangle analysis. At such high precision, 
choices of renormalization scheme and associated pertur- 
bative higher-order effects become an important source of 
uncertainty. Two standard methods have emerged: the 
use of momentum-space subtraction schemes that can 
be nonperturbatively implemented on a lattice [5] and 
the Schrddinger- functional method [3], where so-called 
renormalization-group-invariant (RGI) masses and ma- 
trix elements are obtained via a direct implementation of 
the renormalization group on the lattice. Within the for- 
mer approach, parameters need a further conversion to 
purely perturbative schemes such as MS [I] , where short- 
distance QCD and new-physics effects are best tractable. 

It has recently been realized that the standard 
RlW/MOM prescription suffers from a strong sensitiv- 
ity to IR effects [5], which has become the dominant 
source of uncertainty on the lattice. This is paralleled 
by unusually large higher-order terms in the perturbative 
conversion factors [3] . A modified scheme with much bet- 
ter IR behaviour has been recently proposed and called 
RI/SMOM [7 . In this work, we study the renormaliza- 
tion of the pseudoscalar (non-singlet) density, which by 
virtue of chiral symmetries is related to the renormal- 



ization of the quark mass, and obtain the next-to-next- 
to-leading-order (NNLO, two-loop) conversion factor al- 
lowing to obtain MS light quark masses from their coun- 
terparts renormalized in the RI/SMOM scheme, or its 
variant RI/SM0M 7(j , as 'measured' on the lattice. We 
find much smaller perturbative corrections than in the 
RlW/MOM case, extending one-loop findings in [7] and 
implying percent-level uncertainties on the MS masses. 



RI'/MOM, RI/SMOM, AND RI/SMOM^ 

In the RI'/MOM renormalization scheme for the quark 
field and mass, two conditions [2] 
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are imposed on the inverse quark propagator S^ 1 = 
Z~ 1 S B ~ 1 . The bare quark propagator Sb is defined 
through (our notation closely follows [7]) 

-iS B (p)= Jd 4 xe^(T^ B (x)i; B (0))), (3) 

and the traces are over colour and Dirac indices. 
and |2| determine the renormalization constants Z q and 
Z m relating bare and renormalized field and mass, ipn = 

1/2 

Zq tp B and m B = Z m mB ■ Both renormalization con- 
stants depend implicitly on the regulator (lattice, dimen- 
sional regularization, etc.) and on the gauge coupling 
and the gauge parameter. A virtue of the RI'/MOM 
scheme is that it can be implemented nonperturbatively 
on the lattice as well as in dimensionally regularized con- 
tinuum perturbation theory. The RI'/MOM field and 



2 



mass can then be converted perturbatively to the MS 
scheme via V>F = (zf^/zf /mom )1/3 ^u'/mom and 

m^ s = Z^ s /Zm^ MOM m^ MOM , where all renormal- 
ization constants have to be computed with the same 
(but otherwise arbitrary) regulator. Both conversion fac- 
tors are known to three-loop accuracy [U[S]. However, 
the perturbation series does not converge well, and this 
constitutes a drawback of using the RI'/MOM scheme for 
extracting light-quark masses from lattice simulations. 
Another issue is the influence of non-perturbative long- 
distance physics. This is most clearly seen by considering 
(non-singlet) axial-current Ward identities such as 

5'i3 1 (p')75 + isSg 1 ^) + i(m U! B+m SjB )Ap t B(p,p'), (4) 

where q = p — p' , and the bare vertex functions B for 
the axial current and App for the pseudoscalar density 
are defined through 

Sb(p')A% b (p,p')S b (p) (5) 
d 4 x d 4 y ' x e-™ (T([iu B ^ 7s s B }(0) u B (x) s B (y))) , 

S B (p')K P . B (p. / p')S B {p) 

= ( d 4 x d 4 y e^' x e~^ (T([iu Bl5 s B }(0) u B (x) s B (y))) . (6) 



Q holds for a regulator which respects chiral symmetry 
(in the limit m B — > 0) . This is the case for certain lattice 
regularizations and for dimensional regularization with 
anticommuting 75. (The use of anticommuting 75 is un- 
problematic here as ^ and the formulae below do not in- 
volve closed traces containing odd powers of 75.) To pre- 
serve Q under renormalization, the axial current must 
not be renormalized, and the renormalization constant 
Zp of the pseudoscalar density must satisfy Zp = Z" 1 , 
where Zp can be fixed by imposing the condition 

a«(pW) = z-^a^pVV) 

= Z- 1 Z P tr[Ap iB (p,p')7is]=!=12 (7) 

at a suitable subtraction point. The choice p 2 = p' 2 — 
— /i 2 , q 2 = corresponds to |2]). But at g 2 = 0, 
App(p,j/) receives contributions from the kaon (pseudo- 
Goldstone) pole, which diverge in the chiral limit rap —> 
[2], and is sensitive to condensate effects suppressed 
only by (Aqcd/a*) 2 [!]. In [7], a modified renormaliza- 
tion scheme, termed RI/SMOM, was proposed, which is 
less sensitive to these effects. In that scheme, ^ is im- 
posed at the symmetric point p 2 = p' 2 = q 2 = — fi 2 . 
Following [7J, we will consider a more general kinematic 
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define conversion factors 
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where a B (p 2 ) = ti[S B (p)$\- The right-most expres- 
sion in ([9]) has a straightforward perturbation expansion. 
Moreover, in [7J a variant scheme RI/SM0M 7ti was in- 
troduced where the field-renormalization condition ([T]) is 
replaced by the requirement 



A«(pVV) = Z-i\ B (p 2 ,p'\q 2 ) 
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Z- 1 tr[A^ B (p,p')757A l ]=48, (10) 



which implies conversion factors 
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The schemes for field and mass are converted as 
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where X = RI/SMOM or RI/SM0M 7(J . 

We note that C x and C x depend on ln/x 2 /^ 2 = lnr, 
where v is the dimensional renormalization scale, and 
implicitly on v through the scale dependence of a s and 
the gauge parameter £. Setting /i = v allows relating 
the anomalous dimensions in the RI/SMOM schemes to 
those in the MS scheme [HHH] according to 
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Here we use the definitions (which conform to [7j) 
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with Y = MS or RI/SMOM or RI/SMOM^. 

NNLO COMPUTATION 

We now compute the conversion factors to 0(a 2 ) in 
dimensional regularization (d = 4 — 2e). Let us denote 

cr = -4N c p 2 + +(T (2) +O(o2), (18) 
A = 4A C + A^ 1 ' + A^ + 0(al), (19) 
A = 4rfA c + A( 1 )+A^+C(a^), (20) 

where the superscripts denote the loop order, cr*- 1 -', A^, 
and A^ 1 ' have been evaluated in [7j. For the present com- 
putation, we also need their 0(e) parts, which will affect 



the 0(a 2 ) results for C q and C m . Taking the traces and 
employing partial fractions, we obtain 



5(1,1) 



a 

47T 
1-S 



5(2,1) 



a s d-{l-0 



xiq 2 j(l,l,l;p 2 ,p' 2 ,q 2 )+g(l,l)e^ 



(21) 
(22) 



\%\p 2 ,p'\q 2 )^2N c C F 



Pis 
4 7T 



(23) 



{[(d-2)V-(d-2)(i-0(p 2 +p' 2 )]j(i,i,i;p 2 ,p' 2 ,9 2 ) 



-2(rf-2)(l-0 5 (l,l)e' yEe 



,2 \ e 



2(1-05(1,2) 



+ ((d-2) 2 + (d-4)(l-0)5(l,l) 



-P 



-P 



where je is the Euler-Mascheroni constant, v the dimen- 
sional renormalization scale and 
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The function j results from a massless triangle, via 
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with j<3 = — (pi + p 2 )- Several cases have been evaluated 
in |13j (our j is essentially their J) , in particular 
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where x = p\jp\ and y = p\jp\- The functions (x, y) 
and have been given in [13] in terms of poly- 

logarithms up to second and third order, respectively. 

At the two-loop level, the relevant diagrams are shown 
in Figure [T] They can be represented in terms of three 
master "topologies" (Figure [2]), which may be called 
"propagator" , "ladder" , and "non-planar" , with their 
propagators raised to general integer powers. For the 
latter two topologies, irreducible numerators occur. The 
set can be reduced by standard reduction techniques and 
a systematic application integration-by-parts (IBP) iden- 
tities. For this we employ the program FIRE [15] , a pub- 
lic implementation of Laporta's algorithm |16j and the 
method of S-bases |T7j. A subtle apsect of the IBP re- 
duction is the occurence of quadratic and simple poles 
in e in the coefficients of the resulting integrals. In a 





FIG. 1. Two-loop propagator and vertex diagrams. The grey 
blobs indicate a sum over all one-loop corrections to a prop- 
agator, the black boxes an insertion of a fermion bilinear 





FIG. 2. Basic three- and two-point topologies: ladder, non- 
planar, propagator (from left to right) 



two-loop computation, this leads to poles of up to fourth 
order. On the other hand, the Feynman diagrams have 
poles of at most second order, entirely of ultraviolet ori- 
gin. The spurious third- and fourth-order poles cancel, 
which constitutes a check of the computation, but they 
also imply a possible dependence on terms up to C(e 3 ) 
in the e expansion of the master integrals remaining after 
the reduction. In practice, we find that only known mas- 
ter integrals [TH [H] are needed, except for the unknown 
0(e 2 ) part of j(l, 1, 2 + e). Denoting 
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-ln.Tlny-3(l-x-y)$ (1) (.x,y) + e 2 (x, y) + 0(e 3 )) , 
we find that 

S« (x, x) = ^ (2) (x, x) - Q {2) (l, ^ - (3-6.t)*« (x, x) 
+ yln 3 .T + 14C(3) + lnx^(3-6a;)$ (1) (a;,a;)-^7r 2 ^ , (27) 
S« = -~rt 2 \x,x)+3^ 1 \x,x) 

(28) 
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The function il^ arises in evaluating ladder master in- rithms. Combining all terms and MS-renormalizing the 
tegrals [131 114] and is given there in terms of polyloga- gauge coupling and gauge parameter, we obtain 
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where r — fp/is 2 , n/ is the number of quark flavours, 
and <f>( 2 ) is given in terms of polylogarithms in |13j . The 
function $W has dropped out of the final results. We do 
not know the origin of this cancellation, involving many 
different terms, including the 0(e) one- loop terms. As 
an elk test, setting fj, = v in (29 1 and taking u j — > 0, we 

recover C™' /MOM [8] . The 0(a s ) terms in Q-Q 

agree with [7] (for w = 1 and r = i ). T he most general 
^ri/smom 7m 



form C„ 



from (30) and (31) as 
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The mass and field anomalous dimensions in the two 
schemes ar e ea sily o btai ned by substituting the expres- 
sions H, H, and Jill, as well as C SM0M = c„ RI ' /MOM 



[51 [8] and the well-known two-loop /?- function into (14) 
and (jl5[). More explicitly, denoting 
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0(0) = (11JV C - 2n/)/12, frV = (34 N 2 - lQN c n f - 
6 C pUf ) /48, and the remaining coefficients can be read 
off p^HpTll we have to NNLO: 
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which vanishes in the Landau gauge, is again straightfor- 
ward to evaluate to O(a^) from ( 29 )-( 31 1 and the pertur- 
bation expansion of S defined in (16). Gauge invariancc 



implies 5 = 7a$, where 7^ is the anomalous dimension 



of the gluon field (defined analogously to (16)) 
giving 
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PHENOMENOLOGY 



To explore the phenomenological consequences of our 
result for QCD with three dynamical light quarks (as in 
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FIG. 3. Conversion factors Cm I ' SMOM as function of cj 



(32) 1 If at LO (top/red), NLO (middle/blue), and NNLO (bot- 



tom/black), and c£ I,/MOM = C7* I/SMOM (0) at NNNLO (dot). 




FIG. 4. Residual matching-scale dependence of the conversion 
factor (^m/ smom at a; = 1 at LO (red), NLO (blue), and 
NNLO (black). 



nature, and in modern unquenched simulations), we set 
rif = 3. Figure [3] shows the conversion factor C m (uj) 
in the Landau gauge. We observe that the NNLO cor- 
rection, like the NLO term, is very small at the SMOM 
point oj = 1. This is in contrast to the RI'/MOM scheme 
u) = 0, where even the next-to-next-to-next-to-leading- 
order (NNNLO) correction JS] [5] is large (dot in the Fig- 
ure). To estimate the effects from uncomputed 0{a z s ) 
terms, we vary the renormalization scale (matching scale) 



v used in the conversion and evolve G 



RI/SMOM 



(w = l;i/) 



to the fixed scale /i = 2 GeV, which gives a formally v- 
independent number |T5, 20J. The result is shown in Fig- 
ure^ The width of each band, due to the uncertainty on 
a s (M z ) = 0.1184±0.0007 [18], is almost negligible. This 
is a consequence of the smallness of the NLO and NNLO 
corrections. We observe that the NNLO result is almost 
scale-independent. Alternatively, we can convert the MS 
mass to the RGI quark mass employing the relevant ex- 
pressions in [20], which is also scale-independent. The 
result is similarly stable under scale variation, but the a s 
dependence is a bit more pronounced. A slightly larger 
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residual scale dependence is found for the RI/SM0M 7 
scheme. Numerically, we obtain 
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corresponding to a perturbative uncertainty of less than 
2%, or about 2 MeV for the strange quark mass, when 
converting from the RI/SMOM scheme, and about 3% 
for the RI/SM0M 7 scheme. As also the absolute 
size of the NLO and NNLO corrections is larger for 
the RI/SM0M 7(j scheme, we advocate the use of the 
RI/SMOM scheme together with an appropriate error 
estimate in extracting results for the light quark masses. 



CONCLUSION 

We have computed the RI/SMOM -> MS and 
RI/SM0M 7 — > MS conversion factors for the quark 
mass to NNLO and shown that the RI/SMOM and 
RI/SM0M 7(1 schemes, designed to reduce sensitivity to 
low-energy non-perturbative physics, are perturbatively 
very well behaved, too. These schemes thus may be used 
to extract quark masses with percent-level accuracy from 
numerical lattice QCD. An important question is whether 
the same holds true for other quantities of interest, such 
as Bk and other hadronic matrix elements. 
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and Andrzej Buras for comments on the manuscript. M. 
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Note: After the initial submission of this manuscript 
to the Arxiv, Ref. [24] appeared, whose authors com- 
pute the conversion factor G m at NNLO for the sym- 
metric renormalization point w = 1, where they confirm 
our result. They also give the corresponding field and 

„ „RI/SMOM 7 „ .JU/SMOM-, 

mass conversion factors C q and C m 

for the RI/SM0M 7jj scheme, as well as expressions for 
the NNLO anomalous dimensions in both schemes. In 
this revised version, we have given expressions for those 
quantities, as well. Specialising to uj = uj' = r = 1, our 
results for the agree with the results in [21]. Setting 



further £ = 0, we agree with the results for 7^ given 
there, up to a global sign difference [25] . 
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